Local density of states induced by anisotropic impurity scattering in a d-wave 

superconductor 
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We study a single impurity effect on the local density of states in a d-wave superconductor 
accounting for the momentum-dependent impurity potential. We show that the anisotropy of the 
scattering potential can alter significantly the spatial dependence of the quasiparticle density of 
states in the vicinity of the impurity. 
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A search for the symmetry of the cuprate superconduct- 
ing state involves among others theoretical studies of the 
disordered systems 0, E3, 0] . The most direct probe of 
a simple defect impact on superconductivity is provided 
by the scanning tunneling microscopy (STM) measure- 
ment of the position-dependent quasiparticle density of 
states. The STM images of Zn and Ni substitutions at 
the planar Cu sites in Bi2Sr2CaCu20s+,5 reveal a dis- 
tinct four-fold symmetry of the local density of states 
(LDOS) 0, 0| predicted for the d-wave superconductor 
response to disorder 0,0. In addition to providing 
detailed information on the superconducting state, this 
kind of experiment may shed light on the nature of the 
quasiparticle scattering centers. Although the main fea- 
ture of the four-fold symmetry of the tunneling currents 
is captured by the model of isotropic impurity scatter- 
ing the observed spatial dependence of the quasiparticle 
density of states is far more complex, and may originate 
from a non trivial structure of the impurity potential. 
This aspect of the quasiparticle scattering process has 
been raised in the discussion of macroscopic measure- 
ments in these compounds H H El El) . It has been 
shown that the unexpected for the d-wave superconduc- 
tivity weak impurity-induced suppression of the critical 
temperature can be understood within a scenario of a 
momentum-dependent (anisotropic) impurity scattering 
potential 0, 0, E3. Eq| . In the present paper we dis- 
cuss the effect of anisotropic impurity potential on the 
local density of states in the vicinity of a single impurity 
and verify to what extend this scenario can reproduce 
the STM maps of P^S^CaC^Os+a- There are studies 
of the momentum-dependent impurity scattering which 
indicate the existence of the impurity-bound states even 
for the potential strength far from resonant E3- There- 
fore, we focus here on the symmetry of the LDOS and its 
spatial dependence. For this purpose and for analytical 
simplicity we consider scattering in the Born approxi- 
mation. Hereafter we show that the anisotropy of the 
scattering potential alters the local density of states in 
the d-wave superconductor but preserves its tetragonal 
symmetry. Compared to the effect of the isotropic im- 
purity potential it enhances the spatial variation of the 
quasiparticle density of states and can lead to a long- 
range spatial modulation of the LDOS spectra. We note 
that the anisotropy of the impurity potential, if the same 
as the anisotropy of the superconducting order parame- 
ter, enhances the relative magnitude of oscillations in the 



LDOS for the direction of a gap maximum, especially in 
the vicinity of the impurity, whereas the impurity po- 
tential with maxima in the nodal region inverts this ten- 
dency and leads to a larger quasiparticle density of states 
along the gap node direction. 

We consider the momentum-dependent impurity po- 
tential m miii e m 



v(k,k') = [v l + v a f(k)f(k')}r 3 



(1) 



where Vi, v a are isotropic and anisotropic scattering am- 
plitudes, respectively, / (k) is the anisotropy function, 
and fj (i — 1,2,3) are Pauli matrices. We assume 
/ (k) = ±1 such as its Fermi surface (FS) average van- 



ishes, i.e., (/) 



dS k n(k)f(k) = 0, where n(k) 



FS 



is the normalized angle resolved FS density of states, 
dSkn(k) = 1. Potential determined by / (k) = 

FS 

sgn (k%. — ky) = sgn (cos2(j)), which is in phase with the 
d-wave superconducting order parameter leads to a par- 
ticularly moderate suppression of the critical tempera- 
ture [12, El, E3 • I n the above is a polar angle. 

We study the effect of the impurity potential ||TJ on 
the quasiparticle states using the one particle Green's 
function of the d x 2_ y 2--w&ve superconducting state in the 



particle-hole notation 

Go (k, uj n ) = [iuj n T 
where A(k) = A (kl 



£ k f 3 -A(k)fi] 



(2) 



3) 



Acos2(f> represents the su- 
perconducting order parameter, io n is the Matsubara fre- 
quency, £k = £k — £f is the quasiparticle energy in the 
normal state, ep is the Fermi energy, and fo is the iden- 
tity matrix. For simplicity we assume a parabolic dis- 
persion and a two-dimensional (planar) superconductiv- 
ity. Neglecting the anisotropy of the Fermi surface we 
concentrate on the feature of coupled anisotropics of the 
order parameter and the impurity potential. The impu- 
rity effect is studied at zero temperature by the analytic 
continuation of (J2J 



G (k, u) = G (k, u) n ) ■ _ 



(3) 



In order to discuss the real space distribution of the quasi- 
particle states we take a two-dimensional Fourier trans- 
form of the Green's function ©, which for the parabolic 
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band £k = k 2 /2m and for positive to smaller than the Fermi energy reads 



3C 
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cf A; 
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Go (k,«) 



2(2tt) 



wto + A (k) fi / ^ 



A 2 (k) 



— ik_r 



T3 I e + — e 



(4) 



where similarly to the three-dimensional case 

1/2 



and 



m 



e F ± ,/w 2 - A 2 (k) 



Im^Jw 2 — A 2 (k) > 0. For negative w one must 
multiply the imaginary part of (@J) by sgn(uj). Here 

a = i. 

We consider the impurity effect in the Born approxi- 
mation, that is, proceed with the perturbative expansion 
of the Green's function up to the second order term in 
the scattering potential v. The linear correction results 
from the impurity-induced change of the Fermi energy 



5G' (k, k', u) = Go (k, u) v (k, k') G (k', u) (5) 

while the quadratic term represents the proper Born cor- 
rection and takes the finite quasiparticle life-time into 
account 



5G" (k, k', w) = Go (k, u) ± B (k, k') Go (k', w) 
where the self-energy is determined by 



t B (k, k0 = I> (k, k") Go (k", La) v (k", k' 

k 



(6) 



(7) 



Both corrections are evaluated for a single impurity 
which means that the self-energy is not determined self- 
consistently. We have neglected changes to the absolute 
value of the order parameter in the vicinity of the impu- 
rity which will not change the symmetry of LDOS images 
but their absolute magnitudes. The total correction to 
the Green's function reads 



SG (k, k',Lo) = 8& (k, k', u) + 5G" (k, k', u) (8) 
and its real space transform is given by 



SG (r, to) = e l ( k ~ k ') r SG (k, k', u) 



(9) 



k,k' 



Function © determines the position dependent 
impurity-induced change of the quasiparticle density of 
states near a single impurity 



SN (r, u) = 7m {6G U (r, u)\ 

7T 



(10) 



We have evaluated the LDOS for equal impurity scatter- 
ing strength in isotropic and anisotropic channel irNoVi = 
■nN^Va =0.1, where Nq is the density of states per spin 
at the Fermi level in the normal state, and for the co- 
herence length £o = 12-Kkp 1 . The distance scale in the 
figures is set by kp 1 . We observe a pronounced spa- 
tial variation of the quasiparticle density of states for 
anisotropic impurity potential. This feature is in agree- 
ment with the result of an enhanced LDOS due to a finite 
range of the scattering potential [l^. In Figs, la-c we 
show the distance dependence of the impurity-induced 
change to the LDOS at the quasiparticle energy u> = 1.1 A 
along the a-axis (0 = 0), i.e., gap maximum, and nodal 
(cf) = it/ 4) direction for isotropic scattering (Fig. la) 
and two representative anisotropic impurity potentials 
JTJ: / (k) = sgn{k 2 — k 2S ) — sgn(cos2(j)) which is in 
phase with the d-wave order parameter (Fig. lb), and 
/ (k) = sgn (k x k y ) = sgn (sin2(f>) (Fig. lc) scattering 
out of the order parameter phase. We note that the im- 
purity potential in phase with the order parameter en- 
hances the magnitude of oscillations in the LDOS for the 
direction of a gap maximum, especially in the vicinity of 
the impurity, whereas the out of phase impurity potential 
inverts this tendency and leads to a larger quasiparticle 
density of states in the nodal regions. Significant is also 
very weak anisotropy of the LDOS near the isotropic im- 
purity (Fig. la) where discernible differences between 
4> = and 4> = 7r/4 directions occur at a distance of 
about SOfc^ 1 from impurity It shows that visible spa- 
tial distribution of intensity maxima and minima in the 
STM maps is induced by the anisotropy of defect poten- 
tial. Therefore, we suggest that the apparent anisotropy 
of the experimental images 0, || may possibly result in 
part from the presence of the momentum-dependent scat- 
tering potential. We have also performed the LDOS cal- 
culation for quasiparticles below the gap threshold energy 
and displayed them for to = 0.1A in Figs. 2a-c. Another 
interesting feature of the anisotropic scattering is a pres- 
ence of a long-range spatial modulation of the LDOS for 
the impurity potential in phase with the d-wave order pa- 
rameter, i.e., / (k) = sgn (k 2 — A: 2 ). Such a modulation 
is absent for isotropic or out of phase scattering (Fig. 3). 
We summarize our discussion by showing in Figs. 4-5 
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comprehensive pictures of the STM images at the fre- 
quencies u> = 1.1A and lo = 0.1A around simple defects 
of isotropic and anisotropic Born potential in the ab plane 
of the d-wave superconductor. The potential strength, 
coherence length and distance units have been fixed as 
in Figs. 1-3. The quasiparticle density of states is given 
in the units of the FS two-spin density of states N(0), 
N(0) — 2Nq, and varies from the lowest (black) to the 
highest (white) value according to the scale in each figure. 
We have checked that similarly to the effect on the criti- 
cal temperature |13|| any other scattering potential which 
is orthogonal to the order parameter(in the sense of the 
FS integral as the scalar product) gives qualitatively the 
same LDOS as the discussed above / (k) = sgn (sin.20), 
i.e., out of phase potential. The effect of potentials in 
phase with the order parameter like / (k) = cos2(f> agrees 
qualitatively with / (k) = sgn(cos2<f>). Therefore, both 
studied anisotropies can be considered representative for 
the impurity potential QJ. Although we have restricted 
our study to nonmagnetic impurity scattering, as a clos- 



ing remark we note that inclusion of spin S — 1/2, 1 
or 3/2 scattering does not contribute any quantitative 
change to the results for the impurity potential obeying 
the Born approximation. 

Concluding, we have shown that the anisotropy of the 
impurity potential enhances the quasiparticle LDOS in 
the vicinity of the impurity and leads to discernible dif- 
ferences of its intensity in horizontal and diagonal direc- 
tion. Depending on its symmetry the impurity potential 
can rotate the LDOS maxima and minima by 45° de- 
grees and can induce a long-range spatial modulation of 
the quasiparticle density of states but cannot change the 
induced by the d-wave superconductivity four-fold sym- 
metry of the STM images. 

Note added in proof: A similar modulation to the one 
in Fig. 3(b) is seen along the gap nodal direction for the 
potential out of the order parameter phase. 
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FIG. 1: Distance dependence of the impurity-induced change 
to the LDOS, 8N(r), at the quasiparticle energy u> = 1.1A 
along the a-axis (0 = 0- solid line) and diagonal direction 
(<f> = 7r/4 - dashed line) away from the impurity located at 
r = for a) isotropic impurity potential; b) potential in phase 
with the order parameter, /(k) = sgn(k^. — ky); c) potential 
out of the order parameter phase, /(k) = sgn(k x k y ). JV(0) 
is the FS density of states in the normal state and &f is the 
Fermi momentum. 
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FIG. 2: Distance dependence of the impurity-induced change 
to the LDOS at the quasiparticle energy lu = 0.1 A along the 
a-axis (0 = 0- solid line) and diagonal direction (<f> = tv/4 - 
dashed line) away from the impurity located at r = for a) 
isotropic impurity potential; b) potential in phase with the 
order parameter, /(k) — sgn(k\ — ky); c) potential out of the 
order parameter phase, /(k) = sgn(k x k y ). 



6 

0.010 | , , , 1 — . 1 — . — , — , — , — . 1 



SN (r) 
N(0) 




20 40 60 80 100 120 140 

r M 

(a) 



0.0-5 



SN (r) 




r M 

(b) 

FIG. 3: Long-range distance dependence of the impurity- 
induced change to the quasiparticle density of states along 
= direction at the frequency uo — 1.1 A for a) isotropic 
impurity potential and the out of the order parameter phase 
potential, /(k) = sgn(k x k y ); b) potential in phase with the 
order parameter, /(k) — sgn(k^ — ky). 




FIG. 4: LDOS-map at the frequency ui — 1.1A around the 
Born impurity located at (0,0) in the ab plane of the d-wave 
superconductor for a) isotropic impurity potential; b) poten- 
tial in phase with the order parameter, /(k) = sgn(k x — k y ); c) 
potential out of the order parameter phase, /(k) = sgn(k x k y ). 
The density of states is given in the units of N(0) by the scale 
next to each map and the distance is measured in kp 1 units. 
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FIG. 5: LDOS-map at the frequency ui = 0.1A around the 
Born impurity located at (0,0) in the ab plane of the d-wave 
superconductor for a) isotropic impurity potential; b) poten- 
tial in phase with the order parameter, /(k) = sgn(k^,—ky); c) 
potential out of the order parameter phase, /(k) = sgn(k x k y ). 
Units as in Fig. 4. 



